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Abstract 

Classical scale invariance is one of the possible solutions to explain the 
origin of the electroweak scale. The simplest extension is the classically scale- 
invariant standard model augmented by a multiplet of gauge singlet real scalar. 
In the previous study it was shown that the properties of the Higgs poten¬ 
tial deviate substantially, which can be observed in the International Linear 
Collider. On the other hand, since the multiplet does not acquire vacuum 
expectation value, the singlet components are stable and can be dark matter. 
In this letter we study the detectability of the real singlet scalar bosons in the 
experiment of the direct detection of dark matter. It is shown that a part of 
this model has already been excluded and the rest of the parameter space is 
within the reach of the future experiment. 


1 Introduction 


Higgs boson was discovered in 2012 at the CERN Large Hadron Collider DM- Since 
then, its properties, such as spin, parity and couplings to the standard model fermions 
and gauge bosons, have been measured and it turned out that they are consistent 
with the standard model prediction. In spite of the success of the standard model 
up to now, however, it is commonly believed that the standard model is not the 
ultimate theory of particle physics. In fact there are lots of unsolved problems in 
the field of particle physics as well as cosmology. 

One of them is the origin of the spontaneous symmetry breakdown of the elec- 
troweak gauge group. In the standard model, the electroweak symmetry is broken 
by Higgs field that has an ad hoc tachyonic mass term. One explanation for the 
tachyonic mass is supersymmetry. In supersymmetric extension of the standard 
model, the negative mass term is induced radiatively. On the other hand, radia¬ 
tive symmetry breaking is possible in non-supersymmetric theory, which is known as 
Coleman-Weinberg (CW) mechanism [3j. In the CW mechanism spontaneous sym¬ 
metry breaking is induced at quantum level from classically scale-invariant scalar 
potential. Although it turned out that the CW mechanism with a Higgs does not 
work for the electroweak symmetry breaking, simple extensions of the Higgs sector 
are known to be phenomenologically viable (see, e.g ., (IH2HI)- 

Recently Higgs properties were studied in a classically scale-invariant standard 
model augmented by an electroweak singlet scalars that form a multiplet of global 
O(N) symmetry [3D]. It was shown that the Higgs self-couplings deviate significantly 
from the standard model prediction. Such feature can be observed as a prominent 
signal of this model at the next-generation lepton collider experiment, such as the 
International Linear Collider (ILC) [3TH33]. On the other hand, it was also shown 
that the singlet field does not get a vacuum expectation value (VEV). Then, other 
Higgs properties are unaffected since there is no mixing between the singlet and 
Higgs. Another important consequence is the stability of the singlet field due to 
unbroken O(N) symmetry. If the reheating temperature of the universe is higher 
than the mass of the singlet, the singlet field is thermalized. Then non-vanishing 
thermal relic of the singlet remains, which can play a role of dark matter. 

In this letter we study direct detection of the real singlet dark matter with O(N) 
symmetry. It was pointed out in a similar framework where the thermal relic abun¬ 
dance of the singlet dark matter is too small to explain the present energy density of 
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dark matter by taking into account the 125 GeV Higgs mmm- This is due to an 
enhanced annihilation cross section caused by a large singlet-Higgs coupling. On the 
contrary, however, the large singlet-Higgs coupling may result in a large scattering 
cross section of the singlet with nucleon. According to the study of Ref. [30], the 
couplings of the singlet with the standard model particles are fixed for the successful 
electroweak symmetry breaking via the CW mechanism, which makes it possible to 
determine the relic abundance and the scattering cross section of the singlet with 
nucleon at a high precision. The scattering cross section of singlet scalar dark matter 
is also discussed in several literature mentioned above. We revise the calculation of 
the spin-independent cross section of singlet scalar particle by adopting the formal¬ 
ism given in Ref. [34] where next-to-leading order QCD effect is properly taken into 
account. It will be shown that part of the model has already been excluded by recent 
LUX result [33] and the future experiments will be able to probe almost the entire 
parameter space of the model. 

Here is the organization of this letter. In Sec. [2] we briefly explain the model, 
including the prescription how to determine model parameters. Then the thermal 
relic and the scattering cross section of the singlet are calculated, and the detection 
of the singlet scalar is discussed in Sec. [3] Sec. Q] is dedicated to conclusion. 

2 The Model 

In the framework with classical scale invariance it is known that the standard model 
without the Higgs mass term has already been excluded. In order to construct 
phenomenologically viable model, therefore, it is necessary to extend the model, 
e.g., by adding a new particle to the model. The simplest extension is to introduce a 
gauge singlet real scalar held. Such a singlet scalar can couple to the Higgs in general, 
then the singlet contributes to the CW potential. The effect strongly depends on the 
degree of freedom of the singlet held. To see the impact we introduce a fundamental 
representation of a global O(N) symmetry, S = (Si, ■ ■ ■ , Sn) t . Consequently, the 
tree-level scalar potential which is allowed under the symmetry is 

U = A h (H ] H) 2 + A nsH'HSiSi + 5 A) 2 , (2.1) 

where H is the Higgs doublet held H = (H + , H°) T , and summed over % — 1, • ■ ■ , N 
for N > 2. Z 2 symmetry is assumed for IV = 1 case, whereas it is also a subgroup 
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of O(N) symmetry and always survives in the scale-invariant tree-level potential for 
N> 2. 

The electroweak symmetry breaking is induced via the CW mechanism. To see 
this, the scalar fields can be taken without loss of generality as H = (l/\/2) (0, 0) T 
and S = (ip, 0, • • • , 0) T where 0 and (p are classical fields of the real scalars. Then 
the effective potential at one-loop level is given by 


keff tp') ktree(0; ’p) T I4-loop(0) p) , (2.2) 


with 


T r (, \ , 4 . ^HS ,2 2 , 4 

kfcree(0, P) = —</> + P + ~^P , 

V-Ioop (<t>, v) = jp-jl n < M iV, f) In ^ ^ “ c ‘ 


(2.3) 

(2.4) 


in MS-scheme with renormalization scale fi. Index i denotes the helds which run 
in the loop diagrams, (n;, M 2 and q are given in Appendix |A]) The electroweak 
symmetry is spontaneously broken if dV e s/d(j)[^u) = 0 with (</>) 7^ 0, which implies 

A h ~ - jj^pz as a necessary condition. Then A h should be regarded as the 

next-to-leading order in terms of the order counting of the dimensionless couplings. 
Consequently we rewrite the effective potential as 


Aeff — Vlq + 14 


7 NLO 


(2.5) 


with Vlo and Vnlo being regarded as leading order (LO) and next-to-leading order 
(NLO) of the scalar potential; 

A, 


rr AtfS ,2 2 , A s 4 
Vlo = P + -jP , 
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where F + app is given in Appendix [A] In Ref. [30] it is shown that the successful 
electroweak symmetry breaking with the Higgs mass rrih — 125 GeV can be real¬ 


ized in a given number of N. Table Q] shows the resultscj Roughly speaking, the 


Higgs mass is expected to be m h ~ (\/ 'N\hs/^)vh with vh = 246 GeV, which is 
consistent with the numerical results in Table [T] With the proper order counting, 
the effective potential around the VEV is obtained by replacing the scalar fields as 
4> —> vh + h , </? 2 —> SiSi and expanding by powers of h and s^s,; 


V eS = const. + i m h 2 h 2 + ^ m s 2 + - 


2^2 i x _ 2 „ „ i /v hhh _ U 3 , A hhhh , 4 


■ + 


4! 


A 


1 '7 , 

+ —— Vh h SiSi + 


Xhhss ,2 1 

——h SjSj + 


ssss / \2 , 

SjSd + 


4! 


( 2 . 8 ) 


where m s is the mass of singlet. We have taken /i — Vh and omitted irrelevant 
terms in our later discussion. The results for X hhhh , X hhh , X hhss , X hss , X ssss and m s 
are summarized in Table [3 [3D] Basically \ H and A hs are chosen to give rise to 
vh = 246 GeV and mh = 125 GeV, which determine the couplings (except for the 
singlet self-coupling) and the singlet mass. As, on the other hand, has little impact 
on these results. Since the Higgs self-couplings Xhhhh and A hhh significantly deviate 
from the SM prediction, the precise measurement of the Higgs self-couplings is a 
viable way to test this model. 


Another important fact shown in Ref. [30] is that the singlet does not get VEV 


#3 


Without a VEV of the singlet, Higgs properties, such as Higgs production or decay 
rates, are unaffected. On the other hand, unbroken 0(N) symmetry forbids s* to 
decay. Such stable particles can change the thermal history of the universe. If the 
reheating temperature in the early universe is higher than the singlet mass, the singlet 
particles are thermalized and their number densities freeze out eventually. Then the 
thermal relics can be components of dark matter. In this model the parameters 
which determine the interaction of s* with the standard model particles (the Higgs 
held in our case) are completely fixed as discussed above. Therefore its nature is 


# lr rhis is the results derived from the potential referred as (I) in Ref. [30] . 

# 2 The couplings are obtained from the parameters shown in Tabled] which corresponds to case (I) 
in Ref. (30j. Since the couplings change by a few % in cases (II) or (III), we will use the couplings 
from case (I) in our later calculation. 

# 3 This fact is guaranteed to all orders in perturbative expansion. Strictly speaking, non- 
perturbative effect might break O(N), which could allow non-zero VEV for the singlet. Possible 
(or known) non-perturbative effect is anomaly. In our model, however, O(N) multiplet is scalar, 
thus it is anomaly free. Though one may concern another unknown non-perturbative effect, the 
situation is the same for the standard model, i.e. unbroken t/(l) em symmetry. 
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N 

1 

4 

12 

R 

v H 

= 246 GeV 

Vt 


0.919 


9 


0.644 


9' 


0.359 


A H 

-0.11 

-0.0045 

0.075 

A HS 

4.8 

2.4 

1.4 

A s 

0.10 

0.10 

0.10 


Table 1: The input parameters of the analysis in Ref. [3JTj. The values Ag = 0.10 are 
benchmarks, but there are very few Ag-dependences for A hhh, ^hhhh, A hss- 


highly predictable, and in fact we will see in the next section that the experiments 
of direct detection of dark matter provide the powerful tool to probe the model. 

Before closing this section it is worth noting that N — 1 case is favored by two 
reasons [SD]. In terms of Veltman’s condition [3Sj the level of fine-tuning for the 
Higgs mass at the electroweak scale gets milder compared to larger N cases. Second, 
the fine-tuning is relaxed compared to the standard model and larger N case in a 
sense that the cutoff scale due to Landau pole is predicted to be around TeV. We 
will come back to this point later. 


3 Detection of Singlet Scalars 


As we discussed the singlet scalars are stable and they can play a role of dark matter. 
Since the singlet scalars interact with Higgs boson, they are thermalized in the early 


universe if the reheating temperature is higher than the singlet mass 


#-i 


Then the 


relic abundance is determined by the conventional freeze-out scenario. What we 
need is the annihilation cross section of s*. Relevant annihilation processes are 
SiSi —$■ W + W~, ZZ , hh and tt. The cross sections for the processes are given by 


VsiSi^-W+W- 


(3 f (s,m w 2 ) ( A hss ™w 2 

4v Tsfa V s- m h 2 


2 + i/ i -W 

4 \ m w z 


2 


(3.1) 


# 4 To be specific, the reheating temperature should be larger than 300-600 GeV (see Table E), 
which is a canonical case in the early universe. 
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N 

1 

4 

12 

Xhhh 

1.32 

1.32 

1.32 

^hhhh 

2.9 

2.9 

2.9 

^hss 

11.4 

5.02 

2.80 

^hhss 

14 

5.6 

3.0 

^ssss 

6.5 

1.9 

0.9 

m s [GeV] 

556 

378 

285 


Table 2: The results about the interactions among the Higgs boson and singlet scalar 
bosons derived in Ref. m- They are defined in Eq. (12.811 . Some of them are also input 
parameters of the calculations in Section [3l These values are those in case (I), except for 
A S sss which is in (II) of Table 4 in Ref. [3D]. In our notation of the SM Higgs potential the 
predicted self-couplings are = xj^l = 6A h = 3m/, 2 /r// 2 ~ 0.78. 


OSiSi^tZZ 


&SiSi—>hh 


(3f(s,m z 2 ) (\hs S m z 


8ns Pi 


2 \ 2 


s - m h ' 


1 (s — 2 m z 
2 T — 


2 \ 2 


4 V m z 2 


I6ns(3i 


of 2\\2 , 4AA fes v h . 4 

Pf{s, m h )A +-—-log — 


sfii 


2X hss \ H 4 ( s(3ipf(s,m 2 h 
sfii V M- 


+ 


2m 


2 , 4 

- lo S- 

, - s t_ 




3 Pf 3 (s,m t 2 ) ( A hs S m t 


8n/3i 


S - 77 V 


(3.2) 


(3.3) 

(3.4) 


where mw, m z and m t are the masses of W, Z and t, respectively, s is the center- 
of-mass energy in the initial state, 4 = (1 — 4 m s 2 /s) 1 / 2 , /3f(s,m 2 ) = (1 — 4 m 2 /s) 1 ^ 2 
and 


A = A 


hhss 


+ 


XhhhXhssVH^ 
s — mh 2 


4 = rn h 2 - - [l T PiPf(s,m h 2 )] 


#5 


(3.5) 


It is 


Numerically it is found that the annihilation mode to tt is subdominant 
straightforward to compute the thermal-averaged cross section from above expres¬ 
sion. We use the formula in Ref. 137j to get the relic abundance. (We have checked 
that the density parameter by solving the Boltzmann equation agrees with the ap¬ 
proximated result within a few %.) The results are shown in Table [31 (We note 


# 5 This is expected since the amplitude of fermion pair final state is chirality suppressed, i.e ., the 
cross section is proportional to mt 2 /m Si 4 instead of l/m Si 2 . 
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that when N > 2 all scalars si, • • •, sjv become dark matter.) It has figured out 
that the relic abundance Q iSj . of s l is much smaller than that of dark matter, which 
means that s* cannot be the main component of dark matter. This is due to the 
large annihilation cross section enhanced by the large couplings (mainly Xh ss )- On 
the other hand, however, s* can be detected in the experiment of direct detection of 
dark matter, which we will discuss below. 

For evaluation of the spin-independent cross section of s* with nucleon, we adopt 
the formalism given in Ref. [31]. (See also Refs. [38l I4T] for earlier works.) In the 
present case only scalar-type operators are induced by Higgs-exchange diagram, then 
the effective Lagrangian for the scattering process is 

C, s = Y, CiO's . (3.6) 

i=q,G 

where 


2 - 


O q s = m q sfqq 


Cf = — 

7T ^ 


(3.7) 


m q is quark mass, G^ JU is the gluon field strength and a s is the strong coupling 
constant. By integrating out the Higgs boson (and top quark), the Wilson coefficients 
at the electroweak scale /i\y — fnz at the next-to-leading order in a s are given by 


{nw) 

Cs i^w) 


X 


hss 


2 m h 2 
Xhss 
24m/, 2 



An 


(3.8) 

(3.9) 


The amplitude is given by the hadronic matrix elements, i.e. ( N\m q qq\N ), (N\^Si 2 G' 
{N = p,ri), which are obtained from lattice simulations [421143] and the QCD trace 
anomaly 


#6 


(N\m q qq\N) = m N f^ ] , 

ry Ary ^ 

(N\—G a G afJ,y \N) = m N 


n 


7T/5, 


(TV/=3) 


1 - (1 - 7m) f- 


(N) 


q=u,d,s 


(3.10) 

(3.11) 


Here is nucleon mass, is mass fractions, e.g., ffjP = 0.019(5), = 0.027(6) 

and f'jP = 0.009(22), which are evaluated in Ref. [4T] based on Refs. [421143] . /3 S and 
are the beta function of a s and the anomalous dimension of quark mass defined 

#6 It is given by 0^ = + (1 - yJE q m q qq, and m N = (N\Q^\N) to derive 

Eq.dTITl). 


G atu, \N) 
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N 

1 

4 

12 

dm 

2.01 x 10~ 4 

4.54 x 10 " 4 

8.07 x 10 - 4 

dgj } [ 10~ 46 cm 2 ] 

6.77 

25.6 

74.5 


Table 3: Fractions of the energy density of s, in the observed dark matter density, and 
spin-independent cross sections of the scalars with proton multiplied by their fractions 
in the present dark matter density. The cross sections are computed with the effective 
couplings /^|nlo shown in Table [1J 


by j3 s = ft and 7 m m q = respectively. In Eq. (13. 11 j) the number of flavors 

Nf = 3 is taken since we only know the mass fractions for the light quarks, which 
means that the Wilson coefficients C| and C$ should be evaluated at the hadronic 
scale /ihad — 1 GeV. This can be done by the matching procedure at each quark 
threshold (bottom and charm quarks) and by solving the renormalization group 
equation for the Wilson coefficients. (For the details, such as the matching and 
renormalization group evolution at the next-to-leading order in a s , see Ref. H). 
Finally the spin-independent cross section is given by 


a {N) ~ 
°si — 


m N 


7T (m s + m N ) 2 


l/ S 


W|2 


(3.12) 


with 


4 N> = ^ C|(// hI<i )(iV|m,9 9 |iV)+C s G ( Mhlul )(iV|^G^G“'“'|iV) . (3.13) 

z—/ 7T ^ 

q=u,d,s 


If renormalization group evolution is ignored as well as taking the leading order 
threshold matching, then the effective scattering amplitude is simply given by 


A 


(AO 


A 


hss 


m N 


est 


2 m h - 


2 7 
9 + 9 


E / 

q=u,d,s 


(N) 


(3.14) 


which is often used in the literature. To see the impact of the proper matching 
procedure to the effective coupling f$ N \ we show the numerical values in Table HI 
The difference between the LO and the NLO results is about 4%, while the one 
between the rough estimation and the NLO is about 7%, which gives rise to about 
14% deviation in the spin-independent cross section. 

Now we are ready to see the experimental consequence of the model. Since the 
scalar fields are not the main component of the present DM, the “effective” cross 
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N 

1 

4 

12 

idlest [10 -5 GeV -1 ] 

9.07 

3.99 

2.23 

#|lo [10 -5 GeV -1 ] 

10.2 

4.49 

2.50 

fs^ nlo [10- 5 GeV x ] 

9.77 

4.30 

2.40 


Table 4: Amplitude for Sj-proton scattering. /g^| es t is given by Eg. (I3.14P and /^Ilcu 
/s^Inlo are the results obtained by appropriate matching at the leading order and the 
next-to-leading order, respectively. 


section of the singlet with nucleon is obtained by multiplying the fraction of the total 
abundance of the scalars in the present DM density, 


cr, 


(N) _ X N ) 


si 


= cr, 


SI 


E 

i= 1 ,- ,N 


x 


n 


DM 


(3.15) 


where Ddm = 0.264 [45]. The results are given in Table [3] and Fig. |T| In Table [3] we 
have used the NLO result. Compared with the most stringent bound for the cross 
section (33], N — 12 case has already been excluded at 90% C.L. The others, i.e. 
N = 1 and 4 cases, are still viable, but N = 4 case is very close to the present bound. 
The cross section is far below the present bound for N = 1 case. However, it is 
much larger than the neutrino background in the direct detection experiments [46] . 
Therefore, our model (with any number of N) can be tested in ton-scale future 
experiments, such as LZ program m- Recall that IV = 1 is favored in terms of 
Veltman’s condition as well as the fine-tuning of the Higgs mass, which is discussed 
in the previous section. Thus the result shows that the most well-motivated case 
will be able to be examined in the future experiments. 


4 Conclusion 

In this letter we have studied direct detection of singlet scalar dark matter in a 
classically scale-invariant extension of the standard model. The model extends the 
Higgs sector to have an additional electroweak singlet scalars that form a multiplet 
of global O(N) symmetry, and the electroweak symmetry is broken via Coleman- 
Weinberg mechanism. Recently the Higgs self-couplings as well as new couplings 
and the singlet mass were precisely computed in Ref. [30]. In the work it was shown 
the Higgs self-couplings deviate from the standard model prediction significantly, 
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Figure 1: Spin-independent cross sections of the singlet scalars compared with the LUX 
90% C.L. bound [35]. Orange, green and purple points represent N = 1, 4 and 12 cases, 
respectively. Triangle and circle points represent rough estimation and the next-to-leading 
order calculation, respectively. 


which can be observed at the next-generation collider experiments such as the ILC. 
Another important outcome of their analysis is unbroken O(N). Consequently the 
singlet scalars are cosmologically stable and can play a role of dark matter. Since all 
couplings and mass parameters are fixed for given number of N [30] , it is possible to 
precisely predict the nature of the singlet scalars. Therefore detection of the singlet 
scalars is complementary for the test of the model. 

Assuming that the reheating temperature is above the singlet mass, we have 
computed the thermal relic abundance of the singlet scalars, and the scattering 
amplitude of the scalars with nucleon. For the precise determination of the scattering 
cross section we have used the formalism [3Tj which takes into account the next-to- 
leading order QCD effect in consistent way. We have focused on three benchmarks, 
i.e. N — 1, 4 and 12 (the singlet masses are predicted as 556, 378 and 285 GeV, 
respectively). Then it has been found that although the relic abundance is much 
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smaller than the present dark matter ( 0 Sj /fiDM ~ (9(10 -4 )), the scattering rate is 
enhanced due to the large Higgs-singlet coupling. To be concrete, N = 12 case has 
already been excluded by the LUX experiments, meanwhile N — 4 case is near the 
bound. In N = 1 case which is favored in terms of the fine-tuning regarding the 
Higgs mass, the effective spin-independent cross section (~ 6.8 x 10 -46 cm 2 ) is far 
below the current bound. It is, however, much larger than the neutrino background. 
Thus it is concluded that the whole parameter space of this scenario is testable in 
the future ton-scale detector of dark matter direct detection. 
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A Parameters in the effective potential 

Here is the list of parameters of the one-loop potential in Eq. ( 12 . 21 ) : 

1 5 

nw = 6 , M w ~ = -g 2 0 2 , c w — - ; 

4 6 

n z = 3, M z 2 = - {g 2 + g' 2 )4> 2 , °z = \ ; 

4 6 

n t = -12 , M t 2 = ^Vt 2 (p 2 , c t = \\ 

n± = 1, M± 2 = F±, c± = ^ ; 

3 

^ 5 -1^//\G Y ^HS^P j ^-Hng ^2 ' 

n s NG = N ~ 1 i M Sng 2 = \hs(P 2 + ^sp 2 , c 5ng = - , (A.l) 

where i = W, Z, t show W, Z, t in the loop, respectively, i — ± indicates ip, 0, 
while i = 77nG) ‘S'ng stand for the degrees of freedom which are orthogonal to 0 
and ip, respectively. In the effective potential with the precise order counting (i.e. 
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Eq. (12.5j) ). we use 


E±app(</*) P) 


A 


HS ,2 . A HS + 3Ag 2 

H--- p 


± 



Agg — 3Ag 
2 



+ 4A//g 2 0V 2 > 


(A.2) 


for F ± . The reason for dropping F_ app in Eq. (12.5ft is explained in Ref. [30] (see also 
Ref. [48]). 
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